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Abstract 

We make predictions for the exclusive semileptonic and the non-leptonic decay 
widths of the B c meson. We evaluate the B c semileptonic form factors for different 
decay channels in a relativistic model, and use factorization to obtain the non-leptonic 
decay widths. 



The recent discovery of the B c meson by the CDF Collaboration ffl attracted a great 
deal of attention. The B c meson is very interesting because it carries non-vanishing flavor 
quantum numbers, and lies below the threshold of the BD decay. Therefore, it can only 
decay through weak interactions which makes this doubly heavy meson useful for studying 
the weak decays of heavy flavors. The B c production mechanisms, spectroscopy, and decays 
have been analyzed using different approaches ( see Ref. H for a review ). 

In a previous paper [[J we have used a relativistic model ||] based on the Bethe-Salpeter 
Equation (BSE) to evaluate the spectrum of the B c meson. No free parameters were used 
to fit the B c spectrum. Instead, all the model parameters had been fixed in previous inves- 
tigations of other meson spectra. We also evaluated the decay constant of the B c meson, 
the inclusive decay widths of the c-quark and the 6-quark together with the annihilation 
width. Our results agree very well with the CDF results of the B c mass and lifetime. We 
have presented these results with two covariant reductions of the BSE and observed little 
dependence on the choice of the reduction especially in the heavy flavor sector. 

In this paper we evaluate the exclusive semileptonic B c — ► P(V)eu and two-body non- 
leptonic B c — > PP, PV, VV decay widths, where P (V) denotes a pseudoscalar (vector) 
meson. We use our model to calculate the semileptonic form factors for different decay 
channels. We then use factorization to obtain the non-leptonic decay widths. We will utilize 
primarily a single reduction since this investigation uses BSE results from the heavy flavor 
sector. 

The BSE provides an appealing starting point to describe hadrons as relativistic bound 
states of quarks. The BSE for a bound state may be written in momentum space in the form 

G-\P,p)i;(P, P ) = J -l-Vi^p-p'MP^W, (1) 

where P is the four-momentum of the bound state, p is the relative four-momentum of the 
constituents . The BSE has three elements, the interaction kernel (V) and the propagator 
(G) which we provide as input, and the amplitude {jp) obtained by solving the equation. We 
also solve for the energy, which is contained in the propagator. 

Different approaches have been developed to make the four dimensional BSE more 
tractable and physically appealing. These include the Instantaneous Approximation (IA) 
and Quasi-Potential Equations (QPE) ||. In the IA, the interaction kernel is taken to be 
independent of the relative energy. In QPE, the two particle propagator is modified in a way 
which keeps covariance and reduces the 4-dimensional BSE to a 3-dimensional equation. Of 
course, there is considerable freedom in carrying out this reduction. 
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Earlier, we have used two reductions of the QPE to study the meson spectrum [ffl . These 
reductions correspond to different choices of the two particle propagator used to reduce 
the problem into three dimensions. We refer to these reductions as A and B. Reduction A 
corresponds to a spinor form of the Thompson equation and reduction B corresponds to a 
new QPE introduced in Ref. |7j. These two covariant reductions are chosen because they are 
shown to give good fits to the meson spectrum. In both reductions, we assume the interaction 
kernel to consist of a one gluon exchange interaction, Voge, in the ladder approximation, 
and a phenomenological, long range scalar confinement potential, Vcon given in the form 

Voge + Vcon = 4 a 7m ® 7m I o-lim - *® 1 (2) 

3 (p — p') 2 v->odfi 2 — (p — p') 2 + fi 2 

Here, a s is the strong coupling, which is weighted by the meson color factor of |, and the 
string tension a is the strength of the confining part of the interaction. We adopt a scalar 
Lorentz structure Vcon as discussed in [Q. In our formulation of BSE there are a total of 
seven parameters : four masses, m u =md, m s , m c , the string tension a, and two other 
parameters used to govern the running of the strong coupling constant. We varied these 
parameters to get the best fit for a list of known mesons as described in [|J . 

In our subsequent work |J on the B c meson, we evaluated the B c spectrum without 
changing the parameter values mentioned above ( see Table Q below ) and compared our 
results with those of Eichten and Quigg || and Gershtein at al. @ using both the Martin 
potential and Buchmuller-Tye (BT) potential. The first row of Table [l] should be compared 
with the experimental result [[]J of 6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . We have also 
evaluated the inclusive c-quark and 6-quark decay lifetimes || and obtained a B c lifetime of 
0.46-0.47 ps in good agreement with the experimental B c lifetime of 0.46io;ie (stat.) ± 0.03 
(syst.) ps 0. 

We now turn our attention to exclusive decays. The B c exclusive semileptonic and non- 
leptonic decays have been discussed in the literature || [10, TT[] . The effective Hamiltonian 



for the semileptonic decays has the standard current-current form, and is given by 

Hw = ^^^(1-75)3^(1-75)/- (3) 

The leptonic current is completely known and the matrix element of the vector (V^) and the 
axial vector (A^) hadronic currents between the meson states are represented in terms of 
form factors which are defined by (considering the channel B c — > B S (B*)) 



{B S {P')\V„\B C {P)) = f + (P + p')^ + f_(P-P' 



7<- 



3 



Table 1: Spectrum of B c mesons in different channels (GeV/c 2 ). 



State 


Our work 


Our work 


Eichten and Quigg 


Gershtein et aZ.0 


Gershtein et aZ.0 




Reduction A 


Reduction B 


Ref. § 


Martin potential 


BT potential 




6.356 


6.380 


6.264 


6.253 


6.246 


l 3 ^ 


6.397 


6.415 


6.337 


6.317 


6.337 


l 3 Po 


6.673 


6.692 


6.700 


6.683 


6.700 


1 3 P 2 


6.751 


6.773 


6.747 


6.743 


6.747 


1'PI 


6.752 


6.777 




6.729 


6.736 




6.888 


6.874 


6.856 


6.867 


6.856 




6.910 


6.891 


6.899 


6.902 


6.899 




6.984 


6.955 


7.012 


7.008 


7.012 



{B* S {P' ,e)\V»\B c {P)) = ige^pe^iP + P'riP-P'f, 
(B* s (P',e)\A„\B c (P)) = fel + {e\P)[a + {P + P'), + a_{P-P'),]. 



(4) 



/+,/_, g, f, a-f and a_ are Lorentz invariant form factors which are scalar functions of the 
momentum transfer q 2 = (P — P') 2 where P and P' are the four-momenta of the B c , B s 
(B*) mesons respectively. 

In our formalism, the mesons are taken as bound states of a quark and an anti-quark. 



We construct the meson states as [12 



I-W(Pm) J, mj)) = \/lM J d 3 p(LmLSms\Jmj) (sm s sms\Sms) 



m 



®Lm L (p)\q(—ri- Pm - P,m s ))\q(— — P M + p,m s )) 



M, 



HI 



M, 



(5) 



where the quark states are given by 



\q(p,m s )) 



\ 



(E q + m q 
2m n 



X" 



(E q +rn q ) 



X 



M, 



m a + m, 



E n 



m? + p 2 . 



(6) 



In the above equations M is the meson mass. The meson and the constituent quark 
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states satisfy the normalization conditions. 

(M(P' M ,J',m'j)\M(P M ,J,mj)) = 2E5 3 (P' M - P M )i A; ^, mjI (7) 



E 

(Q(p',rn' s )\q(p,m s )) = — ^ 5 3 (p' - p)5 m > a >ma . (8) 

Tfln 



The wavefunctions $t mi appearing in Eq.(^) for the mesons are calculated by solving 
reductions of Bethe-Salpeter equation H. We have applied this formalism to evaluate the 



semileptonic form factors of the B to D and D* mesons and showed that our results |L3 
are consistent with the heavy quark effective theory (HQET). We use wavefuctions from 
reduction B as we did in our previous work on B decays [|l"3|j . 

The values of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements we use in this 
paper are V ud = 0.974, V us = 0.2196, V ub = 0.0033, V cd = 0.224, V cs = 0.974, V cb = 0.0395 



0- 



In Fig. 1 we show the semileptonic form factors for B c — > B S (B*) and in Table || we 
show the exclusive semileptonic decay widths to different pseudoscalar and vector final states 
(.B+ — > P{V)e + u). We also compare our results with those of |TI[]. We notice in Fig. 1 that, 



although the semileptonic form factors are qualitatively similar to the B — > D(D*) ones JT3 
flavor symmetry is absent in B c decays as discussed in JTH . 



Table 2: Exclusive semileptonic B+ — > P(V)e + v decay widths in 10 6 ev. 





Process 




Decay width 


Decay width 










This work 


Chang and Chen [11 




Bt 


-> i] c e + 


z/ 


11.1 


14.2 


b decay 




* J/ipe 




30.2 


34.4 




Bt~ 


-> 




0.049 


0.094 




Bt~ 


-> D* e" 




0.192 


0.269 




Bt- 






14.3 


26.6 


c decay 


Bt~ 






50.4 


44.0 




Bt' 




V 


1.14 


2.30 




Bt- 


-> 5*°e" 




3.53 


3.32 
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For non-leptonic decays, the effective Hamiltonian (considering the B£ — > B s ir + channel) 
may be written as 

H w = %K.Kd[ci(M)Oi + (9) 

where 

Ox = (uidi) V -A(sjCj)v-A, 

2 = (Uidj) V -A{SjCi)v-A, (10) 

with = 1,2,3) denoting color indices and V — A referring to 7^(1 — 75). C\(fj) and 
c 2 (/^) are short distance Wilson coefficients computed at the scale \x. By factorizing matrix 
elements of the four-quark operator contained in the effective Hamiltonian of Eq.(H), one 



can distinguish three classes of decays [16]. The first class ( class I ) contains those decays 
in which only a charged meson can be generated directly from a color-singlet current, as 
in B+ — ► B s 7i + . A second class of transitions ( class II ) consists of those decays in which 
the meson generated directly from the current is neutral, like the 7r° meson in the decay 
B+ — > B + 7T°. Class I decay amplitudes are proportional to a 1; class II decay amplitudes are 
proportional to a 2 where 

ai = ci(fx) +£c 2 (//), 

a 2 = oiifi) + £c 1 (fi), (11) 

and £ = 1/N C , where N c is the number of quark colors, and \i is the scale at which factor- 
ization is assumed to be relevant. For the third class ( class III ) the a\ and a 2 amplitudes 
interfere. Although the QCD factors a%, and a 2 have been calculated beyond the leading 
logarithmic approximation ||17|| , we will follow the prevailing convention of theoretical pre- 
dictions and express our results in terms of them. As an example the B+ B s ir + amplitude 
takes the form 

A(Bj -> B 3 tt + ) = ^V cs V:M^<^ + \(u i d i ) v ^ A \0><B s \(s j c 3 ) v . A \B c > . (12) 

The matrix elements < B s \{sjCj)v-A\B c > in Eq.(^) have already been evaluated in 
semileptonic decays of the B c meson in terms of form factors, while the other matrix element 
(< 7i + \(uidi)v-A\Q >) is related to the decay constant of the relevant meson. The weak decay 
constants fp and fv for pseudoscalar and vector mesons are defined by 

< 0| J„\P{p) > = ifppp, 

<0\J IX \V( P )> = MyfyS^ (13) 
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where P and V are pseudoscalar and vector states, respectively, and J^, = — is the 
weak current (V^ and are the vector and axial vector currents). The decay constants can 



be expressed in terms of the wavefunctions of the relevant mesons and are given by [ITH 



fi 

F P (p) 
MP) 




p 2 dp 

P 2 



(m q + E q ) (nig + E q ) 

■,2 



1 - 



P 



(m q + E q )(m q + Eg) 

AEgEg 



^v(p), 



Fi(p), 



(14) 
(15) 
(16) 



3(m q + E q ){m q + E q ) 

where ipp(v) are the momentum wavefunctions of the pseudoscalar (vector) mesons. 

We have previously applied this formalism to evaluate the decay constants and the non- 
leptonic decays of the B mesons [19]. The values of the decay constants we use in this paper 
are f n = 0.130 GeV, f p = 0.208 GeV, f K = 0.159 GeV, f K * = 0.214 GeV, f D = 0.209 
GeV, fu* = 0.237 GeV, f Ds = 0.213 GeV, f D . = 0.242 GeV, f Vc = 0.400 GeV, fj^ = 0.400 
GeV. These values are the available experimental ones jjjj]. Otherwise we use our values 
reported in [HI. These values of the decay constants are similar to those used by other 
0,0. 

we compare our results for the exclusive non-leptonic B c —>■ PP, PV, VV 



authors 
In Table 



decay widths of different channels where the b quark decays with those of [11], while in 
Table [|, we make the same comparison for the case of c quark decays. At first glance, our 
decay widths in Table |3| are generally smaller than those of Ref. [O by 20-40%. However, 
this is not a uniform trend as our — > D + D*° is 10% larger than that of Ref. fH|. If 



we furthermore compare total lifetimes for B c we find that our lifetime (0.46 ps) is longer 
compared to Ref. [O] (0.40 ps) which is consistent with the dominant trends seen in the 
comparisons of the exclusive channels. Both theoretical lifetimes are well within current 
experimental uncertainties. Thus, experimental results for a set of exclusive channels could 
resolve between these two sets of theoretical predictions. Table § displays even greater range 



of differences between our model and that of Ref. [II]. 

In conclusion, we have systematically evaluated the decay widths of the exclusive semilep- 
tonic channels B c — > P(V)ev and the exclusive two-body non-leptonic decays B c — > PP, 
PV, and VV assuming that either c or b quark inside the B c meson is a spectator quark, 
and using our relativistic model f4|. In general, our predicted widths are smaller than those 
reported in Ref. |TT| but there are exceptions to this trend. The variations between the the- 
oretical predictions are wide enough so that experimental results should be able to discern 
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Table 3: Exclusive non-leptonic decay widths of the B c meson in 10 -6 ev. b quark decays 
with c quark spectator. The authors of Ref. |ll|] did not report the widths of some of the 
channels because it was thought, prior to the experimental discovery of the B c meson, that 
these channels will be kinematically closed. 



Class 


Process 


Decay width 


Decay width 








This work 


Chang and Chen 11 




C 


—>■ rj c n + 


ajl.59 


a 2 2.07 




C 


— > n r p + 


a 2 3.74 


a 2 5.48 




Bt ~ 


-> J/lfj7T + 

1 T 


a?1.22 


a 2 1.97 




Bt " 


-> J/ibp+ 


a 2 3.48 


a 2 5.95 


I 


Bt 

C 




a?0.119 


a 2 0.161 




Bt- 


-> VcK* + 


a 2 0.200 


a 2 0.286 




Bt~ 


■> J/ipK + 


a 2 0.090 


a 2 0.152 




B c ^ 


J/ipK* + 


a 2 0.197 


a 2 0.324 




Bt - 


-> D+W 


a|0.633 


a 2 0.664 




Bt~ 


-> D+D*° 


a 2 0.762 


a^0.695 




Bt- 

c 


-> D* + W 


a|0.289 


a^0.653 




Bt~ 


* D*+D*° 


alO.854 


a 2 1.080 


II 


Bt ~ 

c 




a 2 0.0415 


a 2 0.0340 




Bt- 


-> DtD*° 


a 2 0.0495 


a 2 0.0354 




Bt- 




a 2 0.0201 


a 2 0.0334 




Bt~ 




0^0.0597 


a 2 ,0.0564 




Bt 




(a!2.16 + a 2 2.57) 2 


(0x1.13 + a 2 1.98) 2 




Bt- 




(a!2.03 + a 2 2.16) 2 


(a!l.04 + a 2 1.90) 2 




Bt~ 




(ail.62 + a 2 1.72) 2 


(ail.02 + a 2 1.95) 2 




B+- 


> J/iJ>D*+ 


(ai3.13 + a 2 3.67) 2 




III 


Bt 




(^0.485 + a 2 0.528) 2 


(ai0.193 + a 2 0.440) 2 




Bt- 


-> r/ cJ D*+ 


(axO.466 + a 2 0.452) 2 


(ai0.181 + a 2 0.430) 2 




Bt~ 


-> J/^L> + 


(ai0.372 + a 2 0.338) 2 


(ai0.177 + a 2 0.442) 2 




Bt~ 




(ai0.686 + a 2 0.732) 2 





between the models. 
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Table 4: Exclusive non-leptonic decay widths of the B c meson in 10 6 ev. c quark decays 
with b quark spectator. 



Class 


Process 


Decay width 


Decay width 








This work 


Chang and Chen [ 


HI 




B + - 


-> B°7l + 


a?15.8 

1 






B+ - 


-> B°o + 


a?39.2 


n^AA 8 








a?12.5 

1 


/-> 2 ^1 (\ 
ti^O 1 .U 




B + - 


-> B*°0 + 


a?171. 


ti-^ 10U. 


I 


B + - 


-> B°K + 


a?1.70 

1 






B + - 


-> B*°K + 


a?1.34 

1 






B + - 


-> B°K* + 


a?1.06 






Bt - 

C 


> BfK* + 

s 


a?11.6 

1 






C 


-> J B°7T + 


a?1.03 


a?3.30 




Bt - 


r 


a?2.81 


a?5.97 




Bt ~ 

C 




a?0.77 


a?2.90 




Bt - 

C 


-> 5* V 


a?9.01 


afll.9 




Bt~ 


^B°K+ 


a?0.105 


afO.255 




Bt ~ 


-> 5°^*+ 


a?0.125 


afO.180 




Bt ~ 


-> 5*°^+ 


a?0.064 


a?0.195 




Bt~ 


> B*°K* + 


a?0.665 


a?0.374 




Bt- 


-> S+^o 


a^39.1 


a^96.5 




Bt~ 


-> B+if* 


a^46.8 


a^68.2 




Bt~ 


-> 


a|24.0 


a§73.3 




Bt~ 


► 5*+^*° 


a^247 


a|l41 


II 


Bt- 


-> £? + 7T° 


a^0.51 


all.65 




Bt- 


-> fi + p 


a|,1.40 


a^2.98 




Bt- 




af.0.38 


ap.45 




Bt~ 


-> £r+p° 


a^4.50 


a^5.96 



We note that the dominant decays are those when the b quark inside the B c meson 
behaves as a spectator quark and a vector meson is produced in the final state. In fact, 
Bt — > B*°e + i> is the dominant decay among all the semileptonic channels ( see Table 2 ) 
and B c — ► B*°p + becomes the dominant among all the two-body non-leptonic decays ( see 
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Table 4 ). Although these decays are suppressed by phase space they are CKM favored. 
Finally we point out that the ratio 

r(B+-+v P +) C _ 2|T , ,™ 



67r 2 \V ud \ 2 aif' (17) 



dr(B+^Ve+v)/dt\ t=m$ 1 - »» 

with V = B*°, J /if> will be a good experimental test for the numerical value of the coefficient 
a x of QCD pj. 
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Figure Caption 

1 The semileptonic form factors for B c — > B S (B*). 
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